Exact Density of States in a System of Disordered Two Dimensional Dirac Fermions 
in the Unitarity Limit: The d-wave Superconductor 



Catherine Pepin and Patrick A. Lee 

Department of Physics, M as sachus setts Institute of Technology, Cambridge, MA 02139 

(February 1, 2008) 

In this paper we present a method to compute the exact density of states induced by N non magnetic 
impurities in a system of two dimensional Dirac fermions in the unitarity limit. We review the case 
of the 7r-fiux phase of the Heisenberg model and also treat the disordered d-wave superconductor. 
In both case we find additional states in the gap with 8p{u>) ~ rii/ \u (in 2 |oj/Ao| + (7r/2) 2 ) |. 



As a consequence of the breakdown of Anderson's the- 
orem [jlj , when impurity scattering violates the symmetry 
of the superconducting state, the superconducting energy 
gap is depleted and impurtites act as strong pair break- 
ers. This is the case in s-wave superconductors with mag- 
netic impurities (^|| for which there is creation of bound 
states in the gap. For the same reasons of symmetry vi- 
olation, non magnetic impurities act as pair breakers in 
unconventional superconductors with higher orbital mo- 
mentum such as d-wave superconductors. As a result, the 
measured low temperature properties of YBaCu207 Q 
displays a remarkable sensitivity to the presence of non 
magnetic impurities: the critical temperature, for exam- 
ple, is suppressed even in the lowest order in the disorder 
potential ||. 

Furthermore the d-wave superconductor is special, due 
to the presence of gap nodes which prevents the complete 
freezing of scattering processes at low energy. 

The standard method to treat these problems of disor- 
der combines the T-matrix approximation with stan- 
dard impurity averaging techniques. For three dimen- 
sional systems such as polar superconductors or heavy 
Fermion superconductors Q the standard perturbative 
approach is reliable. In the limit of low impurity con- 
centration rii, a perturbative expansion in rii leads to a 
finite density of states at the chemical potential |§JlO|| . 

For two dimensional systems (to which it is believed 
the high-Tc cuprates belong) the standard procedure of 
averaging over impurities may be complicated by the ap- 
pearance of logarithmic singularities in the perturbative 
expansion of the single electron self-energy [2C[ |. Such 
a situation appears in a variety of two dimensional sys- 
tems characterized by a Dirac-like canonical spectrum. 
In order to control the diagramatic expansion, all the 
diagrams with crossing impurity lines would have to be 
resummed which is technically impossible. 

In these systems the existence of a finite density of 
states induced by disorder is highly controversial Jl3| . 
Some non perturbative methods have been used to treat 
a weak disorder potential |^,|l],^5). The solution then 
depends on the symmetry of the pure system. In the 
special case of the d-wave superconductor the density of 
states still vanishes at low energy in the presence of disor- 
der. On the other hand some perturbative self-consistent 
calculations |fl 12 and a non perturbative one |l4|] show 



a finite density of states. 

This issue of finite density of states is crucial for the 
conduction properties in the disordered compound. In- 
deed, if there exist states in the gap, possible anomalous 
overlaps between well separated impurities can induce 
a new conduction mechanism entirely through impurity 
wave functions in a so-called 'impurity band' p3| , p4j . 
This arises the question of localisation. 

In this paper, we treat the impurities in the unitary 
limit (infinite scattering potential Vb), but randomly dis- 
tributed in space. In this limit each impurity is identical, 
and create an extended bound state at E = 0. These 
states overlap, but a singular density of states remains 
at zero energy. Our approach is thus very different from 
models of random on-site energies considered by others 
and really does not contradict them. 

We introduce a new method to calculate exactly the 
density of states induced by N non magnetic impurities 
in a system of two dimensional Dirac fermions in the 
unitarity limit. We will outline the general features of 
the proof and present a different derivation than in our 
previous work p2[ . The question of averaging will be 
discussed in details and we show how the result is exact 
provide the average over the different configurations of 
impurities is achieved at the end of the calculation. 

Our result can be applied to several systems of Dirac 
fermions in two dimensions. We review the case of the 
7r-flux phase and explicitely treat the d-wave supercon- 
ductor. 

The generic hamiltonian for a d-wave superconductor 
can be written 



H 



k 



(1) 



with Ek — W (cos k x + cos k y ) — 



and Afe = Ao (cosk x — cosk y ) 

which describes quasiparticles in the presence of the spin 
simglet order parameter Afc. The <7j are the Pauli ma- 
trices in the particle-hole space and ^ = (c\ 
is the spinor corresponding to the creation of a particle 
and a hole. We shall present our results in terms of a 
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d x 2_ y 2 state eventhough our conclusions are general for 
any state where vanishes linearly along a direction 
parallel to the Fermi surface. 

Instead of using the Nambu formalism, we work with 
the diagonalized version of (|l|) in order to access directly 
the properties of quasiparticles. 

The Bogoliubov transformation gives Ckt = Wkak — 
i>k/3k and c ^Lki = ^k^k + u k/3k where and (3^ are 
the annihilation operators for quasiparticle and quasi- 
hole respectively. The hamiltonian then reads H = 

Ek^k («l"k - 0lPk) with uo k = y/e'i + Ai. The co- 
efficients Uk and Vk satisfy u 2 , = 1/2 (1 + Sk/uik), = 
1/2 (1 — Ek/uJk) and Ukfk = Afc/(2cjfc). In order to sim- 
plify the calculation, we consider that the spectrum is lin- 
earized around the four nodes ^±A;_f, ±&f^ where close 

to 1/2 filling, kp — (tt/2 — p) is close to tt/2. 

The disorder is introduced via a scalar scattering 
potential Vo located randomly in the system Hi = 
Vo Ej=i c ! c i; where N is the number of impurities. 

The lagrangian density can be written in the general 
form 

N 

£ kk' = G^ 1 S a p5kV -V ^2 (fLfk'f) + /ka/k'< 

i=l 

(2) 

This describes a system of two dimensional Dirac 
fcrmions in which N non magnetic impurities located 
at Ri, i = 1 to N, have been introduced as repulsive 
scalar potentials of strenght Vo. Quasiparticles with pos- 
itive energy (a — 0) and quasiholcs with negative en- 
ergy ( a = l)are described by the bare Green's function 



T) a Wfc. The factors / kQ and / kQ come 



(, k„ 

from the diagonalization of the initial hamiltonian. Two 
different terms usually arise from this diagonalization, 
leading to the sum in eq. (||). 

In the case of the d-wave superconductor, we have 



ikR; 



and /ka = *kae lk ' Ri where t K o 



fL = (-l) Q tk Q e* 
«k and tk.i = Vk- 

In a previous paper Q we have studied the case of 
the 7r-flux phase, considered as a model for the spin-gap 
phase of the cuprates Jl8| , ^9| , ^2l . Following the same di- 
agonalization as exposed in [p2| , we would get the fac- 
tors: / kce = / kQ = e ik Ri if i belongs to sublattice A and 
fL = fL = (-l) a e ik Rl e i¥,k if i belongs to sublattice B. 

Starting from the lagrangian (^) , the equations of mo- 
tion for the fully dressed Green's function can then be 
written 



^-kq^. 



10 
qk< 



(3) 



We define fka and Gk>/3 as N components vectors in 
the space of impurities of components / ka = Uk a e lk Ri 
and gi, p = E q7 fq-yGqk' and respectively for fL and 

Qk'fl- 



The Dyson equation (without any impurity averaging) 
can be written as: 

G kk' = G l_ a 5kk'8 a 



oGka (fla ' Sk'fJ + PL ' ^Ic'/j) i ( 4 ) 



where the scalar product in the parentheses runs over the 
positions of impurities. 

Solving with respect to Q and Q' gives 



M 



G°f 

g, J \G°f J ' 



(5) 



where M = ( ^ ~t f ^ * 
\ B I + A 

Aij — Eka ^ka/ka/ka 

and B i:j = E kQ G LfLfL( R j) ■ 
We thus get the T-matrix equation 

G kk' = G kc Ak"5a/3 

+ G ia { fka fk a ) T [ ^ j G k> /: 



(6) 



with f = -VoM -1 . 

The additional density of states is given by bp = 

— l/TT-fm-Eka ^kk- Using the equality dG°/diuj n = 
G° 2 and the cyclicity of the trace we get Sp(u>) — 

— 1/irImTr (^M^ 1 dM / diuj^j . The trace is over the po- 
sitions of impurities. 

In our previous study [ p2[ we solved the problem by 
writing the term in the brakets as the derivative of a log- 
arithm and evaluating the resulting determinant. Here 
we use another trick which consists in using the form of 
A/ 2 to factorise the logarithmic divergence caracteristic 
of these problems. This derivation will make clear the 
question of averaging on disorder. We thus write 



5p{uj) = —ImTr \hr 2 dM 2 / 'diu 7 
2ir L 



(7) 



and carefully study the form of M 2 . 

The solution of the problem depends on the actual 
form of My . 

As for the 7r-flux phase, we can show that My w 1 /Rij . 
This fact was previously acknowledged by other oth- 
ors |23|,^5| who found an overall factor of the order of 

(Ao/WW/Rij- 

Let's first consider the case of symmetric bandwidths 
for positive and negative energies. This means we are 
at half filling for the superconductor (p = 0). We show 
below how this l/R value arises in the case where W = 
Ao- 

We get 
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Ai 



E - 



2 c ikR; 



keBZ 



lUJ n — UJk 



E 

keBZ 



so that 

= E 



keBZ 



E 

kesz 



iuJn + UJk 



k II. 



lUJ n + UJk 



= 2nT°(R lJ )-^\n 



A 



2in 



F (Rij) , 



(8) 



with F°(R) = 2 [cos f + i? y ) + cos f (i^ - Ry)] and 

^ rl (i?) = 2 sin -^(Rx + R y ) (cos ip + sin ip) 
+ sin -^{Rx — Ry) (cos ip — sin ip) 

where ip is the angle between R. and the x axis. 
For B we write 



B 



i j 



E 

keBZ 

so that we find 



u k i; k e 



ik R 



u k v k e 



k-R, 



kesz 



iw„ + UJk 



B, 



2l7T .-, . „ . 



(9) 



with 



F 2 (R) = 2 sin - (R x + R y ) (cos p - sin p) 

7T 

+ sin — (i? x — i? a ) (cos 93 + sin ip) 

The key point in the evaluation of the density of states 
is that M 2 shows a logarithmic divergence. 

We factorize explicitely this divergence by taking 



M = In 



An 



s 



(10) 



where S is a matrix which depends on the particular 
configuration of impurities we work with. Two terms 
emerge in the density of states: 



Sp(uj) = ImTr 

7T 



hr 



A 



ainK/Ao| ^ dS 



duj 



diui n 
(11) 



Note that the first term doesn't depend on the particu- 
lar configuration of impurities. The average over disorder 
will act only on the second term. After analytic contin- 
uation, wc find per unit of volume 



5p(uj) 



Hi 



uj (In 2 |w/Ao| + (tt/2) 2 ) t:V 



4-Jm TrS-^ 



where the brakets denote the average over disorder. 

We argue now that the contribution coming from S is 
negligible. For this purpose, we make an overestimate of 
S with the N/2 x N/2 matrix S a satisfying S?. = 1 if 
Rij < Ao/|w n | and elsewhere. In this approximation 
we have overevaluated the long distance intercations be- 
tween impurities by replacing the slow decay of which 
vanishes at R^ = |An/iu> n | by a step function with a cut- 
off atRij — \D/iu> n \. This special form of S thus leads 
to an upper bound for the averaged second term in ( |l2| 
). Our aim is now to prove that it is negligible as com- 
pared to the first term. dS a /dui is then the matrix which 
elements dS a / duj\ij vanish everywhere except on the ex- 
ternal boundary of the circle of radius iZy = |A /w n |. 

The main difficulty in the inversion of S a is that two 
circles of radius A /|a>„| centered around two points i 
and j very close to each other will overlap, leading to the 
same number of non zero coefficients in the lines i and j of 
S a . In order to differentiate the sums J2k ^ik^ii' 1 an( l 
Ylk^ik^kj~ we nave thus used the external boundary 
of the circle to compensate its volume in S a _1 . 

We take £?■ ~* = (\w n \/(y/nA )) U tj . For all j inside 
a circle of radius Ao/|w n | around the point i, Uij is a ran- 
dom configuration of ±1 so that Ylj Uij ~ \/7rAn/|a; n |. 
In addition all the points j situated in the external 
boundary of this circle have Uij = —1/(2^/tt). Elsewhere 
= 0. 

We thus get S a _1 (dS a /duSj = ^/n wich is negligible 

as compared to the first term. Thus this term is of second 
order as compared to the first term in (|l2|). 

Provided M 2 displays a logarithmic divergence, we can 
factorize this divergence and show that the exact density 
of state in the unitarity limit and in the limit of low 
frequencies is 



6p(uj) 



|^| (In 2 | W /Ao| + (tt/2)s 



(13) 



In the case where Ao ^ W the logarithmic singularity 
gets an overall prefactor A /W p3[ which leads to the 
most general form at low energy 



5p(uj) 



An 



W 



(14) 



M (ln 2 | W /A | + (V2) 2 ) ' 

The equation ( |l4| ) is thus the most general expression 
for the density of states induced by non-magnetic impu- 
rities in the unitarity limit, and for energy bands sym- 
metric for positive and negative energies. 

Note that this expression is normalisable. If the im- 
purities were totally uncorrelated, we would find N 5- 
functions at uj = 0. Here we see clearly the overlap of 
impurity states which leads to a broadening of the 6 func- 
tions. The interacting states repeal each other, but not 
enough to kill the divergence. This level repulsion was 
found by other methods pQ ] for Dirac fermions in two di- 
mensions in presence of disorder, but away from unitarity 
limit. 
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What happens if the bandwiths are not symmetric any- 
more, that is if fi ^ but still /i< Ao ? 

First the nodes are moved away from the point 
(±7r/2, ±7r/2) so that transversal nodes are now sepa- 
rated y the vectors Q = (ir(l — <5),7r(l — 5)) and Q* = 
(— 7r(l — <5),7r(l — 5)) where 5 = /i/Ao and /i is the in- 
crease in the chemical potential. This leads to a change 
of the phase factors in Aij and B^j, Namely we get 
J^{R) = 2 [cos (Q ■ R/2) + cos (Q* ■ R/2)]; 

T 1 (R) = 2 [sin (Q • R/2) (cos ip + sin p) 

— sin (Q* • R/2) (cos tp — sin <p)] , 

and 

T 2 {R) = 2[sin(Q-R/2)(cos(^-sin^) 

— sin (Q* • R/2) (cos tp + sin <p)] . 

As 6 is a small parameter, this change in the phase 
won't affect the existence of the logarithmic divergence 
in M 2 . 

More importantly, away from half filling, the bands 
of quasi particles and quasi holes become asymmetric 
to account for the removing of particles in the system. 
This gives additional terms in and B y of the form 

their 



J2i<\k\<i+s ellc ' R /( iuj n - 0Jk)- We show that 
contribution in is of order of 5Jo(Rij)/ Aq where Jo 
is the Bessel function of rank zero. This doesn't affect 
the logarithmic divergence at low energies. 

The density of states is thus stable with respect to a 
small perturbation of the chemical potential. 

The result (14) holds for any system of Dirac fermions 
in two dimensions, provided the hopping matrix between 
inpurities goes like 1/Rij. Thus it seems to be generic 
for two dimensional Dirac fermions in presence of Poisso- 
nian disorder. Indeed the same result as ( |l4| ) has been ob- 
tained for a system of electrons in a strong magnetic field 
in two dimensions, when the disorder is Poissonian [ p7| . 

This result has to be compared with one dimensional 
spin-Pcicrls and two-leg ladder compounds for which 
non magnetic impurities create a Dyson-like singularity 



5p{w) ~ 1/ (M In 3 \lu/D\) In one dimension, the 

Dyson singularity is a signature of delocalised states in 
the system: there is a direct relation between the localisa- 
tion lenght and the density of states |28j which shows the 
divergence of the localisation lenght at low energies. In 
two dimensions this theorem doesn't hold anymore and 
the question of localisation is still open. 
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